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ABSTRACT 

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission- 
line galaxies at z ~ 2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extra- 
galactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with > 5cr 
detections of emission lines to / > 2.5 x 10” 18 erg s” 1 cm” 2 , means that the galaxies in the sample 
are typically ~7 times less massive (median M* = 1O 9 5 M 0 ) than previously studied z ~ 2 emission- 
line galaxies. Despite their lower mass, the galaxies have [Oill]/H/3 ratios which are very similar to 
previously studied z ~ 2 galaxies and much higher than the typical emission-line ratios of local galax- 
ies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the 
two-dimensional spectra of the galaxies for this purpose. In the stacked data the [O ill] emission line is 
more spatially concentrated than the H/3 emission line with 98.1% confidence. We additionally stack 
the X-ray data (all sources are individually undetected), and find that the average T[om]/To. 5 -iokeV 
ratio is intermediate between typical z ~ 0 obscured active galaxies and star-forming galaxies. To- 
gether the compactness of the stacked [O ill] spatial profile and the stacked X-ray data suggest that 
at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei. 

Subject headings: galaxies: abundances — galaxies: active — quasars: emission lines — galaxies: 
evolution 


1. INTRODUCTION 


Based on observations with the NASA/ESA Hubble Space 
Telescope , obtained at the Space Telescope Science Institute, 
which is operated by AURA Inc, under NASA contract NAS 
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Rest-frame optical spectra contain a wealth of informa- 
tion about galaxy properties. In particular, the strengths 
and ratios between collisionally excited and recombi- 
nation emission line s can be used to infer star forma- 
tion rate (SFR, e.g. lKennicuttlll998lh gas-phase metal- 
licity (e.g.lKewlev et al.l 200 111, and nuclear activity (e.g . 
I Baldwin. Phillips fe Terlevichlll98lt IKewlev et al .1120061) . 
However applying these techniques to galaxies at higher 
redshifts is challenging for two reasons: (1) at z > 1.5 
optical emission lines are observed in the near-infrared 
(NIR), where ground-based spectroscopy is challenging, 
and (2) the characteristics of z > 1 galaxies may be phys- 
ically different from local galaxies. 

The Wide Field Camera 3 (WFC3) G141 slitless 
grism on the Hubble Space Telescope ( HST) addresses 
the first problem, with far improved capabilities for 
NIR spectrosco py than ground-based observatories (e.g., 
iStrauehn et al.l 1201 If) . Most of the previous ground- 
based studies of z > 1.5 galaxies have b een biased 
to massive and UV- luminous systems jShaj3le£ i et_al 


2004: lErb~ et al. 2006; Havas hTet al.ll2009l : lOnodera et al 


20101 ). probably missing most of the more numerous, less 
massive galaxies at 1.5 < z < 3. Only a handful of 
lower-mass galaxies at z ~ 2 have been studied, revealing 
much lower metallicities than both z ~ 0 galaxies of the 
same mass and th eir massive and evolved counterparts at 
the same redshift (lYuan & Kewlevll2009HErb et al.ll2010l 
iFinkelstein et al.ll2011fl . In contrast the tremendous sen- 
sitivity of HST/ WFC3 allows for > 5cr detections of emis- 
sion lines to / > 2.5 x 10” 18 erg s” 1 cm” 2 for every ob- 
ject in a 2'.l x 2'.l field with 8 orbits. In addition, HST is 
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WFC3 G141 Grism 



Fig. 1. — Line luminosity limits for detecting the Hct, H/3, [O III] , 
and [Oil] emission lines with 8 orbits of the HST/WFC3 G141 
grism (corresponding to a line flux limit of 2.5 X 10 -18 erg/s/cm 2 ). 
Dashed vertical lines indicate the redshift ranges at which each line 
is observed in the G141 wavelength range (1.1 < A < 1.7/rm). In 
this paper we study galaxies with H/3 and [O III] lines, selected at 
1.3 < 2 < 2.4. 

not limited by atmospheric emission or absorption in the 
NIR, and so spectral lines are detectable over the entire 
1.1-1. 7 /i m wavelength range of the WFC3 G141 grism. 
Figure [T| shows the detectability of emission lines with 
the G141 grism at various redshift ranges. 

Star formation processes were qui te different at z > 1 
than in the loca l universe (e.g. IPapovich et al.l 120051 : 
iReddv et al.ll2006l : lShim et a, 1.11201 ll) . and so it is not clear 
if the emission-line diagnostics used to characterize local 
galaxies can be directly applied to z > 1 systems. For 
example, the bulk of z ~ 2 galaxies have higher ratios 
of collisionally excited to recombin atio n li nes ( [O ml /H/3 
and [Nil] /Hot) than local galaxies dShaplev et al.l[2005bt 
lErb et ahl 120061 : iKriek et al.l l2007fl . Local galaxies with 
high SFR tend to have higher emission-line ratios than 
galaxies with lower SFR, suggesting that rapi d star for- 
mation leads to different ionization propert ies (ILiu et al.l 
120081: IBrinchmann. Pettini fe Charlodl2008D . Since z > 1 
galaxies typically have high er SF R than local analogs of 
the same mass (e.g. IPapovich et al.l 1200511 , their higher 
line ratios may simply be the res ult of their higher 
SFR. However [Wright et al.l (|2010T ) used spatially re- 
solved spectroscopy of one galaxy at z = 1.6 to show 
that these line ratios increase in the center of the galaxy, 
with normal star- forming (SF) line ratios at outer aper- 
tures. This suggests that weak or obscured active 
galactic nucleus (AGN) activity may be the cause of 
higher line ratios in z ~ 2 galaxies, perhaps because 
AGN activity is correlated with SFR a s observed locally 
(IKauffmann et al.ll2003HShi et al.ll2009ll . There are many 
predictions for the density of obscured and Compton- 
thick (CT, Nh > 10 24 c m 2 ) AGNs at 1.5 < z < 3 from 
the X -ray background dGilli et al.1 l2007t iTreister et al.l 
120091) and from IR-exces s galaxy counts and X-ray stack- 
ing (IDad diet al. 1 120071 IM artinez-Sansigre et al.l I2007t 
iFiore et al.l 120091 IGeorgantopoulos et al.l 120111) . but to 
date only three C T AGNs have been ind i vidually con- 
firmed at z ~ 2 ([Alexander et all 120051 : lErlund et al.1 
120081 [Alexander et alJl2008D . Obscured AGNs are missed 
by even the deepest X-ray surveys, but can still ex- 


hibit th e narrow emi ssion-line signature of an AGN (e.g., 
lLaMassa et al.ll20TIl) . Measuring the density of obscured 
AGNs at z ~ 2 is especially important to reveal how 
much black hole growth occurs in an enshrouded phase. 

Here we use 8.5 orbits of WFC3 G141 spectroscopy to 
study the properties of emission- line galaxies at 1.3 < 
z < 2.4 in the Hubble Ultra Deep Field (HUDF, 
iBeckwith et al.l [200611 . These data were taken as part 
of the supernova follow-up program in the Cosmic As- 
sembly Near-Infra red Deep Extragalactic Legacy Sur- 
vey ( CANDELSP^l. IGrogin et al.ll20lil : IKoekemoer et al.1 
120111) . The supernova c andidate identi f ied in the HUDF 
region is described by iRodnev et al.l (I2011D . and this 
work makes use of the spectra that were serendipitously 
obtained for hundreds of other targets in the 2.1' x 2.1' 
field. In particular we focus on the 28 galaxies exhibit- 
ing H/3 and [O ill] emission lines in the observed grism 
waveleng th ra nge (a t 1.3 < z < 2.4). Like the study by 
IWright et al.1 (I2010D . the high spatial resolution of the 
HST/W FC3 similarly allows for studies of spatial gradi- 
ents in emission-line properties. Unlike an integral field 
unit spectrograph, however, the WFC3 grism has the 
unique ability to do this for hundreds of objects in the 
field of view simultaneously. 

In §2 we describe the observations, data reduction, and 
emission line luminosity limits. In §3 we describe our 
method for calculating emission lines, complicated by the 
limited spectral resolution of the WFC3 G141 grism. In 
§4 we stack the spatially resolved spectra to show that 
both the emission-line properties and the stacked X-ray 
data are suggestive of weak AGN activity in at least 
some of these galaxies. We summarize our results in 
§5. Throughout the paper we adopt a cosmology with 
h = 0.70, =0.3, D A = 0.7. 

2. HST/WFC3 SLITLESS GRISM OBSERVATIONS 

The HUDF was observed with the HST/ WFC3 G141 
grism on 26-27 Oct 2010 and 1 Nov 2010 for a total of 8.5 
orbits. The data were acquired in two different orients, 
with 2.5 orbits at an orient of PA = 176° on the sky and 6 
orbits at an orient of PA = 182.4°. The d ata for each ori- 
ent w ere reduced using the aXe software dKiimmel et al.1 
200!). available at http://axe.stsci.edu/axe/), which 
flux- and wavelength-calibrated the data and provided 
co-added two-dimensional (2D) spectra as well as opti- 
mally extracted one-dimensional (ID) spectra. The us- 
able wavelength range of each reduced spectrum is 1.1 < 
A < 1.7/iin, with a resolution of R ~ 130 (46.5A/pixel 
for a point source and 0'.' 13/pixel). 

The G141 grism is slitless, and this introduces a num- 
ber of complexities when interpreting the data. Spatially 
extended sources have features which appear extended 
in the wavelength direction, and so broad emission or 
absorption lines can be caused by either large velocity 
widths or by the size of the object. In addition, spec- 
tra frequently suffer from contamination by 0th order, 
1st order, and 2nd order spectra of nearby sources. For 
this reason we did not combine the two orients, instead 
analyzing each orient independently and using them to 
identify contamination. 

We calculated redshifts for all 455 targets detected 
in the HUDF G141 data. The best redshift was found 

18 http://candels.ucoIick.org 
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using a cross-correlation template-fitting IDL algorithm 
based on the publicly available idlspec2d package writ- 
ten by D. SchlegePd. Previous spectro s copic redshifts 
ll Wolf et, ~aTl 1200 4 Mignoli et al.1 120051 : iVanzella et, al.1 


120081: S traughn et al.l I2008D and photometric redshifts 
dWolf et al.l 12004 iGrazian et, al.l [20001: iRvan et al.l [20071 : 
iWuvts et al.ll200fl iCardamone et al.ll2010fl were used as 
priors. We used two templates: a star-forming emission- 
line galaxy and a passive absorption line galaxy. Each 
template was constructed using the high-resolution op- 
tical (3800 < A < 9200A) templates from the SDSS li- 
brarj0 combined with supplemental low-resolution UV 
and NIR data from the SWIRE photometric templates 
(jPolletta et al.ll2007h , with the final templates resampled 
to match the spectral resolution of the WFC3 G141 data. 
All spectra and their redshift fits were visually inspected 
to ensure quality. 

Of the 455 targets, we selected the 28 galaxies at 1.3 < 
z < 2.4 and H/3 or [O ill] emission-line luminosities of 
L > 1.8 x 10 41 erg s -1 . This luminosity limit represents 
the flux limit for detecting emission lines at the 5<r level 
(/ > 2.5 x 10 -18 erg s _1 ) at z = 2.4. For all 28 galaxies 
we detect both H/3 and [O ill] in both grisrn pointings, 
and as a result their redshifts are unambiguous. Table 
[I] shows the properties of these 28 emission- line galaxies. 
Their median redshift is 1.86. 


2.1. Stellar Masses and Rest-frame Colors 

Rest-frame magnitudes and stellar masses were cal- 
culated from the GOODS-MUSIC photometry catalog 
(|Santini et al.l [20091 by fitting a set of template spec- 
tral energy distributions dra wn from t he PEG A SE s tellar 
population models (see iFioc fe Rocca-Volmerangelll997l 
for a description of an earlier version of this stellar pop- 
ulation mode l). We allow dust attenuation following 
iGalzettil ( 2001 ). with values of E(B — V) gas between - 
0.05 and 1.5. (Slightly negative E(B — V) avoids biasing 
the mean properties of the population by accounting for 
the rare cases when photometric error gives fits that pre- 
fer slightly negative attenuation. In practice all of our 
emission-line galaxies have best-fit E{B — V) > 0 any- 
way.) All templates have solar metallicity, and included 
a broad range of exponentially-decreasing, constant, or 
exponentially- rising star formation beginning at Zf ~ 4. 
A grid of stellar population templates were compared 
with photometric data points of each galaxy, given the 
spectroscopic redshift from the grism data, and the best 
matching template was used to calculate rest-franre col- 
ors and a stellar M/L ratio. Compared to galaxies with 
independently-estimated stellar masses (from indepen- 
dent photometry of simil ar b ut not i den tical datasets) 
in the GOODS-S field ( Wuvts et al.ll200j ). our method 
yields similar masses for intensely star-forming galaxies 
both with and without dust attenuation (i.e. , the galax- 
ies of interest here), with a scatter of 0.3 dex. Our masses 
are 0.1-0.25 dex larger for moderately star-forming galax- 
ies or non-star-forming galaxies because our templates 
include older ste llar populations than those used by 
IWuvts et al . (l200§ ). Rest- frame colors are reproduced to 
within 0.1 mag in rest-frame ( U — B) color for all galaxy 

19 http://spectro.princeton.edu/idlspec2d_install.html 

20 http://www.sdss.org/dr5/algorithms/spectemplates/ 



Fig. 2. — Rest-frame (t/ — B) color and stellar mass for our 28 
emission-line galaxies (black crosses). Typical lcr errors are ~0.1 
magnitudes in ( U — B) and ~0.3 dex in mass. For comparison 
filled gray circles show the 137 galaxies at 1.3 < z < 2.4 in the 
HUDF and gray contours and points show the entire set of 4140 
GOODS-S galaxies in the redshift range, both with masses and 
rest-frame quantities calculated from the same photometry and 
with the same method. Blue dia mond s show a comparison sample 
of z ~ 2 galaxies from Shaolcv ct al. (2005a): this is the parent 
sample of most previous studies of emission-line galaxies at z ~ 2. 
Arrows show the effects on ( U — B ) color if the H/3 and [Om] 
emission lines are removed (although the [On] emission line, which 
is not measured here, could have the opposite effect on ( U — B) 
color). The WFC3 grism reveals galaxies with masses a factor of 
~ 7 lower than previous studies of z ~ 2 emission-line galaxies. 

types. Full details on th e derived masses a nd rest-frame 
quantities are gi ven by iBell et all (lin prepD . 

As shown by Ivan der Wei et all (|201ll ). strong emis- 
sion lines can have significant effects on broad-band col- 
ors in z ~ 2 galaxies. In our galaxies the median to- 
tal equivalent width (EW) of the H/3 and [O ill] lines is 
EWahs = 220A: this brightens the rest-frame B magni- 
tude by ~ 20% and reddens the ( U — B) color by 0.2 
magnitudes. Despite these effects, we do not correct 
rest-frame magnitudes for emission lines for two reasons. 
First, the [Oil] emission line could have the opposite ef- 
fect, brightening U and causing U — B to appear bluer, 
but we do not measure this emission line because it is not 
generally present in our observed spectral range. In ad- 
dition, the SF/AGN diagnostic we use (lYan et al.ll201ll ) 
was calibrated with (U — B) colors uncorrected for emis- 
sion lines, and so we do not correct our colors to remain 
consistent. 

Figure [2] shows the rest-frame (U — B) colors and 
stellar masses for our 28 emission-line galaxies. Also 
shown by gray filled circles are the 137 galaxies in the 
HUDF with photometric or spectroscopic redshifts at 
the same 1.3 < z < 2.4, with rest-frame colors and 
stellar masses computed in an identical fashion. The 
emission-line galaxies studied here are rather common at 
1.3 < z < 2.4, making up ~20% of the total H < 26 
galaxy population. Gray contours additionally show 
the set of all 4140 GOODS-S galaxies with photomet- 
ric or spectroscopic redshifts in the same 1.3 < z < 2.4 
redshift range, and our galaxies have a roughly simi- 
lar distribution of mass and slightly bluer colors than 
the larger parent distribution. For comparison with pre- 
vious sa mples, we also s how th e 50 galaxies at 1.3 < 
z < 2.4 of lShaplev et al.l (|2005al) . using their masses and 
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TABLE 1 

HUDF WFC3/G141 Emission-Line Galaxies 


# 

RA 

Dec 

z 

Hab 

lo §( L H/3) 

lo s( L [Om]) 

f[Om]/fu/3 

(U - B ) rest 

log(AA) 

/o.5— 10keV 

- 

[deg] 

(J2000) 

- 

[mag] 

]log(erg/s)] 

]log(erg/s)] 

- 

[mag] 

[log(Mg)] 

[ 10 -16 erg s _1 cm -2 ] 

1 

53.135250 

-27.781778 

1.44 

22.67 

4i-o8t8:3! 

4i.48l°;i| 

9 CC9+0.77 
z.oz_i is 

0.76 

10.39 

< 0.45 

2 

53.141159 

-27.759295 

1.77 

24.19 

40.891??! 

41.611?;?? 

c: 00 + 1-22 

d.zo_2.06 

0.51 

9.39 

< 0.44 

3 

53.143940 

-27.779781 

1.71 

25.21 

40.9218;?! 

41.351°;°? 

2.971??? 

0.38 

9.05 

< 0.47 

4 

53.144127 

-27.773602 

1.90 

23.58 

41.1818;?! 

41 Q2 + 0 - 05 

5.54!?;?? 

0.36 

9.51 

< 0.44 

5 

53.144485 

-27.791162 

1.42 

23.52 

40.971°;!* 

41 QS - * -0 ' 05 
4J.-yo_o.o5 

8.09!!;?? 

0.46 

9.20 

< 0.46 

6 

53.144840 

-27.778469 

2.22 

24.37 

41.101°;?* 

41 Q 5+ 0 - 02 
4i.yo_o.o6 

7.101!;?! 

1.00 

9.28 

< 0.47 

7 

53.146435 

-27.778353 

1.87 

24.74 

41.181°;!° 

41 94~*”S'Si 
4i.y4 _ 0 09 

5.74l?;°? 

- 0.01 

8.43 

< 0.47 

8 

53.146900 

-27.781837 

1.55 

24.09 

40.741°'?* 

41 27 +0 - 10 

41.Zf_o 13 

q qc+0.97 

O.OQ_i 5i 

0.48 

9.39 

< 0.46 

9 

53.147469 

-27.777693 

1.83 

22.80 

40.831°;** 

41.57!°;°? 

4.52!!;!! 

1.05 

10.43 

< 0.47 

10 

53.147747 

-27.765675 

2.31 

23.75 

41.081°;?* 

4i-7o!8:8! 

4.22!!;°° 

0.59 

9.91 

< 0.43 

11 

53.147804 

-27.771420 

2.17 

24.52 

41.361°;°? 

42 01 + 0 - 04 
4Z.Ul_o.o9 

5.09!?;*? 

0.65 

8.82 

< 0.45 

12 

53.148003 

-27.787737 

1.90 

23.80 

40.961??? 

41.24!°;!° 

o rjQ+0.35 
z.uz_i 09 

0.63 

9.92 

< 0.45 

13 

53.149063 

-27.785116 

2.07 

24.20 

41.021°;?* 

41.38!°;!° 

1 QO "*” 1 ' 28 
l.yU—o .88 

0.65 

9.72 

< 0.44 

14 

53.151283 

-27.792444 

1.86 

24.05 

41.181°;*° 

41.74!°;!? 

Q 62"*" 0 ' 77 

o.uz_i 62 

0.48 

9.59 

< 0.43 

15 

53.151691 

-27.763506 

2.16 

25.19 

40.551°;*° 

4l.24l°;°* 

4-oo_i . 80 

0.48 

8.94 

< 0.44 

16 

53.152264 

-27.770132 

1.85 

23.21 

41.261°;?* 

41 77 + 0 - 06 
‘ 1 -0.04 

q 9 q + 2.34 

O.ZO- 0.88 

0.53 

9.97 

< 0.47 

17 

53.152870 

-27.772545 

1.85 

24.13 

41-291°;°* 

41.95!°;!* 

c co+1-58 
°- oz -1.36 

0.36 

9.16 

< 0.46 

18 

53.152874 

-27.780163 

1.86 

24.33 

40.651?;?? 

41.3218;°! 

q f 7 q + 0.90 

O. to_ 1 72 

0.60 

9.46 

< 0.46 

19 

53.153774 

-27.767345 

2.32 

23.65 

40-981?;?? 

41.5818;!! 

3.94!!;?° 

0.44 

9.77 

< 0.46 

20 

53.154453 

-27.771494 

2.23 

22.76 

41.671?;?? 

42.49!°;?! 

6.63!!;°! 

0.66 

10.24 

< 0.46 

21 

53.155643 

-27.779324 

1.85 

22.10 

41.601?;?? 

41.83!°;!* 

1 9y+1.33 
-0.32 

0.76 

10.72 

< 0.45 

22 

53.156372 

-27.767874 

2.02 

24.07 

40.751?;*? 

41 2S - * -0 ' 12 
41.ZO_o.20 

2 . 6 i!?;?! 

0.83 

9.19 

< 0.47 

23 

53.157207 

-27.778561 

1.31 

23.86 

40.941?;*! 

41 45 + 0 - 08 
41. 4O_0. 28 

3.17!°;°! 

0.54 

9.20 

< 0.45 

24 

53.160435 

-27.775246 

1.41 

20.86 

41.081?;!° 

41.34l°;°I 

i.si!!;?? 

0.42 

9.64 

< 0.46 

25 

53.160446 

-27.790394 

1.61 

23.44 

41.231?;!? 

41.57!°-°? 

1.39!?;?! 

0.48 

9.74 

< 0.43 

26 

53.165157 

-27.781664 

2.22 

24.56 

40.901?;?? 

41.60l8.i9 

5.88!°;?* 

0.69 

9.32 

< 0.44 

27 

53.172153 

-27.763655 

2.21 

24.83 

41.031?:?* 

41.7918;!? 

c 't't+3.09 
' * -1.84 

0.75 

9.04 

< 0.48 

28 

53.172241 

-27.760641 

1.54 

21.83 

41.331?;?? 

41-65!°;!° 

2.09l°;?° 

0.72 

10.75 

< 0.48 


( U—B ) colors calculated from their ph otometry usin g the 
publi cly available kcor rect package (Blanton fc Roweisl 
1200711 . The galaxies of IShaplev et al. (I2005afl form the 
parent samp le for most previous z ~ 2 emission-line 
studies (e.g. IShaplev et akll2005bl : lErb et al]l2006L l2010t 
IWi ight et al.l2010h . The median mass of the grism galax- 
ies is M* ~ 10 9,5 Mm, a factor o f ~ 7 lo wer than the me- 
dian mass of the IShaplev et al.1 (I2005afl sample. Several 
of the emission-line galaxies from the WFC3 grism data 
appear to be star-forming dwarf galaxies with masses of 
only M* ~ 10 9 Mq. 

2.2. X-ray Data 

The HUDF has the deepest X-ray data in the sky, 
with 4 Ms of Chandra data. Imaging was obtained in 
54 observations (obsIDs) over the course of 3 Chan- 
dra observing cycles in 2000, 2007 and 2010 using the 
Advance d CC D Imaging S pectrometer imaging array 
(ACIS-l lGarmire et al.ll2003f) . The data were reduced us- 
ingCIAO vT2 according to the basic procedure described 
in [Laird et all (|2009H . Before combining the observa- 
tions, the astrometry of each obs ID was registered to t hat 
of the GOODS-MUSYC survey (jGawiser et al.l[2006i l by 
matching the positions of bright X-Ray sources to H- 
band selected sources, using the tool reproject_aspect. 
Source detectio n was ca rr ied o ut according to the method 
described in lLaird et al.1 (120091 1 . and none of our galaxies 


are detected. We instead assign 3cr upper limits calcu- 
lated using the background and PSF size from the loca- 
tion of each galaxy on the Chandra / ACIS detector. 

3. LINE MEASUREMENTS 

Because the grism is slitless, any spatial extent along 
the dispersion direction will artificially broaden spec- 
tral features. The [Om]A4959 and [Oiii]A5007 lines are 
often completely blended together in resolved galaxies, 
and larger systems additionally have H/3A4861 partially 
blended with the two [O ill] lines. For this reason we 
do not attempt to individually measure each of the H/3 
[Oill] lines directly from the data. Instead we subtract 
a linear continuum in the region 4770 < A < 5100A and 
fit the residual with a mixture of three Gaussians. The 
three Gaussians are constrained to be centered within 
20l +1200 km s" 1 ) of 486lA, 4959A, and 5007A, in 
the rest frame. Because the two [O ill] lines are blended, 
we measure the [Oiii]A5007 flux as 3/4 of the total [O ill] 
flux (i.e. , 3/4 of the the sum of the two Gaussians cen- 
tered near 4959A and 5007A; iStorev fc ZeippeiJ I2000D . 
Errors in line flux (and /([Oiii])//(H/3) and L[om]) are 
calculated by bootstrapping 10 000 realizations of the 
resampled data. 

Figure [3] shows the H/3 and [O ill] line regions for eight 
emission-line galaxies. The 2D dispersed image, ID spec- 
trum, and line fit is shown for each object. 
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Fig. 3. — Examples of fits to the H/3 and [O ill] emission lines for eight of the emission-line galaxies. Dashed vertical lines indicate the 
line centers for the H/3A4861, [Om]A4959, and [Om]A5007 emission lines. Because the G141 grism has such low spectral resolution, the 
H (3 and [O ill] lines are frequently blended. We subtract the continuum, shown by the red line, and fit the entire system with a mixture of 
three Gaussians, shown by the green line. The [Om]A5007 flux is measured as 3/4 of the total [O ill] flux. 


The classical “BPT” diagnostic 

(jBaldwin. Phillips fe Terlevichl I1981D uses emission-line 
ratios (namely, /([Om])//(H/3) and /([Nil])// (Ha)) 
to distinguish galaxies dominated by star formation 
or black hole activity in the local universe ( see also 
iVeilleux fc Osterbrockl 119871 : iKewlev et al.lfeOOll ). How- 
ever in the WFC3 grism data, Ha and [Nil] are 
hopelessly blended, and additionally are redshifted 
be yond the observed grism redshift range at z > 1.6. 
lYan et aD (|2011l) showed that it is possible to construct 
an alternate version of the BPT diagnostic which 
replaces /([Nn])//(Ha) with rest-frame ( U — B) color. 
This modified BPT diagnostic is especially useful for 
interpreting the WFC3 grism data. 

An alternative SF / AGN classification dia gnostic possi- 
ble for our galaxies is the “MEx” diagram (iJuneau et al.1 
hud, which replaces the (U — B) color with stellar mass 
M*. Stellar mass, however, is dependent on the physics 
of the model, and could be systematically inaccurate due 
to the constraints of available templates. Meanwhile rest- 
frame color does not depend on the physics of the model, 
only its goodness of fit to the data. F or th i s rea so n we 
prefer (U — B) over M*, noting that the lYan et al.l (120111 1 
diagnostic a grees well with the “MEx” method (see Ap- 
pendix A of IJuneau et al.ll2011l) . 

Figure Q] shows our z ~ 2 e mission-lin e galaxies on 
the modified BPT diagnostic of lYan et al.l (1201 ill , along 
with comparison sam ples at z ~ 0 an d z ~ 0.6 (als o 
from lYan et al .1120111) and z ~ 1.2 Hrom lLiu et al.ll2008l) . 
Galaxies dwelling on the left and bottom of each panel 
tend to be SF dominated, while systems in the upper 
right are AGN dominated. Low metallicity galaxies tend 
to have bluer colors and higher /([Om])//(H/3) ratios, 
and indeed the locus of SF galaxies tends to shift to the 
upper left from z ~ 0 to z ~ 0.6 to z > 1. Our galaxies 


generally lie in same region as previously studied z > 1 
galaxies. 

Directly interpreting the position of z > 1 galaxies on 
the modified BPT diagnostic in Figure [4] is challenging 
for severa l reaso ns, however. First, the dividing line of 
lYan et all ([201 11) between SF galaxies and AGN was em- 
pirically calibrated at z < 1 only, and it is unclear if this 
dividing line changes with redshift. Local galax ies wi th 
higher SFR ten d to have higher line ratio s dLiu et al.l 
2008: iBrinchmann. Pettini fe Charlotl I2008D . and since 
z ~ 2 galaxies have higher SFR galaxies th an local galax- 
ies of similar mass (IPaoovich et al.l[2005D . we might al- 
ready expect them to have a higher /([Oiii])//(H/3) ra- 
tios from both their rapid star formation and low metal- 
licity. On the other hand, active galaxies might also be 
more likely to be confuse d as SF-dominated in lo w-mass, 
low-metallicity galaxies. Ilzotov fc Thuanl (|2008i) showed 
that AGN accreting as high as the Eddington limit con- 
tribute less than 10% to the emission-line ratios in metal- 
poor dwarf galaxies. Our z ~ 2 galaxies have similarly 
low masses, and the emission from their low-mass black 
holes would be overwhelmed by a SFR of a few solar 
masses per year. 

For these reason it is difficult to interpret the abso- 
lute position of objects, and more instructive to consider 
relative positions of the inner and outer regions of the 
galaxies, as well as differences in the profiles of [O ill] 
and H/L 

3.1. Stacked HST/WFC3 Data 

The slitless WFC3 grism provides spatially resolved 
spectroscopy, with a scale of 0"13 per pixel. This is 
particularly useful for studying gradients of emission- 
line ratios: most local galaxies have a higher ratio of 
/( [O ill] )/ /(H/3) in outer regions due to negative metal- 
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Fig. 4. — Ratio of [Om] to H/3 with rest- frame (U — B) color. The solid line shows the division between AGN (upper right) and star 
formation (left and bottom) from [Tan et al.| (1201 If). Crosses show the 28 z ~ 2 galaxies from the WFC3 G141 data. Contours and points 
show comparison samples (also from [^m et al.ll20111) a t z ~ 0 (left panel) and z ~ 0.6 (middle panel). In the right panel, gray diamonds 
show a comparison sample at z ~ 1.2 from ILiu et all (120086 . The red points show the colors and line ratios of the stacked data: the 
diamond shows the inner (//39 region and the circle shows the outer C//39 — 0 / /78 region. The z ~ 0 and 2 : ~ 0.6 comparison samples have 
the same line luminosity constraint as the WFC3 G141 data, with either H/3 or [O ill] of L > 1.8 x 10 41 erg s -1 . From z ~ 0.1 to z ~ 0.6 
to 2 : ~ 1.2, galaxies have lower metallicities and conse quently move t o the upper left of the figure, and our z ~ 2 galaxies lie in very similar 
space to the previously studied z ~ 1.2 galaxies of ILiu et all (120086 . The line ratio and color of the inner region of the stacked data is 
marginally (l.Ocr) more AGN- like than the outer region. 


licity gradients (e.g.- lvan Zee et aPll998li . while a higher 
ratio in a galaxy nucleus would suggest the presence of 
an AGN. In general the signal- to- noise (S/N) of our data 
is too low to measure line ratio gradients for individual 
galaxies; instead we stack the 2D spectroscopy to study 
line ratio gradients. 

We stack the 2D images for each galaxy, equal- 
weighting each galaxy by normalizing its total (H/3 + 
[Oiii]A4959 + [Oiii]A5007) line emission. Images are 
aligned by the peak of the emission (dominated by the 
[Oiii]A5007 line). This stacked spectrum allows us to 
study the spatial extent of the H/3 and [O ill] emission 
lines. 

One-dimensional spatial profiles can be constructed by 
extracting each emission line along the cross-dispersion 
direction. Figure [5] shows the spatial profiles for the 
H/3A4861 and [Om]A5007 emission lines, along with the 
collapsed cross-dispersion profile of the entire 2D spec- 
trum. Each emission line is extracted from the cen- 
tral three pixels (140 A) about the line center. The 
H/3 line has a cross-dispersion full-width half-maximum 
(FWHM) of 7.2 pixels (0"94, kpc), while the [Om] 
line is narrower with a FWHM of 4.9 pixels (0'/64, 
~5 kpc). 

We test the significance of the [O III] line being more 
spatially compact than the H/3 line by bootstrapping 
100 000 resampled profiles. We construct each boot- 
strapped profile by adding a normally-distributed ran- 
dom error (with a given by the propagated error of the 
stacked profile) to the original data, and then measure 
FWHMQOiii]) and FWHM(H/3). The fraction of boot- 
strapped profiles with FWHM([Om])jFWHM(H/3) gives 
the significance of [O III] being more compact than H/3. 
This occurs in 98.1% of the resampled data sets, and so 
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Fig. 5. — Spatial profiles of the H/3 (blue line) and [O III] (red 
line) emission lines, along with the entire spectrum (black line, 
offset by +0.4 in flux). Each ID profile is extracted along the 
cross-dispersion direction, using the three pixels (140A) about each 
line center for the emission lines or the entire wavelength range for 
the entire spectrum profile. The H/3 profile is more extended than 
the [O III] profile with 98.1% confidence, suggesting that at least 
some of the objects have a compact AGN contributing to the [Om] 
emission. 

the [Om] line is more compact than the H/3 line with 
98.1% confidence. 

We also measure the spatially resolved /([Oiii])//(H/ 3) 
ratio and (U — B) in the inner and outer regions of the 
stacked spectrum. We designate the “inner region” as 
the central 3x3 pixel box (0"39 x 0'/39, or 3.2 x 3.2 kpc 
at z = 1.86) in the cross-dispersion direction, and the 
“outer region” as the two boxes 3 pixels wide from pixel 
numbers 4 through 6 (0"39 — 0'/78) above and below the 
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Fig. 6. — X-ray luminosity with [O ill] luminosity for the 28 z ~ 2 
galaxies from the WFC3 G141 data. Since all of our galaxies are 
undetected in X-rays, we show them as X-ray upper limits. The red 
diamond shows the stacked X-ray data (and mean Lroill]) f° r the 
sample. The dashed line s hows t he typical relat ion between Lx and 
L[om] for Type 2 AGNs llLaMassa et al.ll20l II) and the dotted line 
s how s the typical relation for star-forming galaxies IIRovilos et all 
200 91 . both at z ~ 0, with the ler scatter about each relation 
indicated by the gray shaded area. The position of the stacked 
X-ray data in Lx — Troiii] space suggests that an AGN is present 
in at least some of the galaxies. 


center in the cross-dispersion direction. Line ratios are 
calculated for both the inner and outer regions, using 
the same method of continuum subtraction and Gaussian 
fitting described above. 

Color gradients were calculat ed using the CANDEL S 
WFC3 Y, J , and H imaging (iKoekemoer et al.1 1201 Il f. 
The WFC3 images are drizzled to a finer resolution 
(0"06) than the WFC3 spectroscopy, and so we adopt 
(y/36 (6 pixel circular diameter) as the inner region and 
0"36 — 0"72 (circular annulus from 6 to 12 pixels) as the 
outer region. Rather than stacking, we measure aperture 
Y, J, and H for each individual galaxy. We then treat 
the median Y — J (for low redshift z < 1.8) or J — H 
(for z > 1.8) difference from the inner (< 0'.'36) to outer 
(0"36 — 0"72) aperture as the change in (U — B) lest over 
the same aperture. 

As Figure |4] shows, the stacked inner region has 
marginally redder colors and slightly higher line ratios 
than the outer region. In other words, the inner region 
is more like an AGN than the outer region. However the 
large errors (particularly on the outer ratio) mean that 
this is significant only at the l.Ocr level (67%). This is 
probably because the 0'/39 size of the inner region corre- 
sponds to the inner 3.2 kpc of the galaxy (at the median 
redshift z = 1.86), and so includes much more than the 
typical ~1 kpc narrow line region (NLR) of an AGN. The 
marginal difference between inner and outer line ratios is 
merely suggestive, but is reinforced by the far more sig- 
nificant (98.1% confidence) result that the [O ill] profile 
is more compact than H/3. 

3.2. Stacked Chandra Data 

In Figure [6] we show the full-band (observed 0.5- 
10 keV) X-ray luminosity upper limits with the 
[Om]A5007 emission line luminosity. The dashed line 


shows the typical Lx — jqom] relation f or [Om]-selected 
narrow-line AGN (lLaMassa et al.ll2011f) , with associated 
la scatter of ~0.5 dex indicated by the gray shaded area, 
and the dotted line shows the typical Lx — Tjjoni] rela- 
tion for star- forming galaxies dRovilos et al.ll2009 j. with 
the gray shaded area similarly indicating the 1 <t scatter 
of ~0.3 dex. 

We stack the X-ray data for the undetected sources to 
learn whether the population has Lx — from] dominated 
by black hole activity or star formation. The stacked 
data is marginally detected (2.9cr significant) with a flux 
of /o. 5 — iokeV = IT x 10~ 17 erg s -1 cm -2 . At the me- 
dian redshift of our sample, z = 1.86, this translates to 
an X-ray luminosity of To. 5 - iokeV = 2.7 x 10 41 erg s 1 . 
Because the stack is only marginally detected in the full 
band, we cannot comment on its hardness. Unless the 
typical Lx — A join] relation for star formation changes 
with redshift, the position of the stacked data is interme- 
diate between the loci of AGN and star formation. This 
suggests that at least some of the galaxies have an AGN 
contributing to the [O ill] and X-ray emission. 

4. NATURE OF Z ~ 2 EMISSION-LINE GALAXIES 

The stacked set of our z ~ 2 galaxies shows a more 
compact [O ill] profile than the H/3 emission line. This is 
the reverse of wh at is expected for typical galaxy metal- 
licity gradients (ivan Zee et al.l [l998h . but is expected if 
an AGN in the galaxy’s center contributes to the line 
emission. We find a marginally higher [O III] /H/3 in 
the center of the stacked d ata, s imilar to the higher- 
confidence result of I Wright et al.l ( 2010 ) , who found a 
steeper [Om]/H/3 ratio in the center of one z = 1.6 
galaxy. In addition the stacked Chandra data suggests 
an AGN component in at least some objects, since the 
position of the stacked data in Lx — Trom] is inter- 
mediate between the lines of star formation and black 
hole activity. For these reasons we suggest that some of 
these z ~ 2 emission-line galaxies, and tho se wit h similar 
line ratios and colors from other authors dShaplev et all 
l2005bt lErb et al.ll2006t iKriek et al.ll2007lf . probably con- 
tain weak AGN in their centers. 

Because we rely on stacked data to reveal an AGN, it is 
difficult to know the frequency of nuclear activity in these 
galaxies. In the local univ erse, AGN act ivity is known to 
correlate with SFR (e.g. iSlii et al.ll2009h . so AGN activ- 
ity might be expected in our galaxies given their bright 
emission lines and presumably high star formation rates. 
While some authors have argued that AGN are rare in lo - 
cal galaxies with M* < 10 10 Mq (jlvauffmann et al1l2003fl . 
recent work suggests instead that A GN incid ence is con- 
stant with host mass at z < 1 (lAird et all 1201 ill . Our 
sample has a median mass of M* = 10 9 - 5 M 0 , and so 
unless a few massive objects are dominating the stacked 
X-ray and grism spectrum, this suggests that AGN are 
common in low- mass star- forming galaxies at z ~ 2. 
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